Transport through a superconductor-interacting normal 
metal junction: a phenomenological description 
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We propose a phenomenological description of electronic transport through a normal 
metal/superconductor interface of arbitrary transparency, which accounts for the presence of 
electron-electron interaction in the normal metal. The effect of interactions is included through 
an energy-dependent transmission probability that is inserted in the expression for the current- 
voltage characteristics of a non-interacting system. This results in a crossover from the Andreev 
to the tunneling limit as a function of the energy at which transport is probed. The proposed 
description reproduces qualitatively the results obtained with formally correct theories as well as 
experimental observations. In view of its simplicity, we expect our approach to be of use for the 
interpretation of future experiments. 
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Low-energy electronic transport through the interface 
between a normal metal (N) and a superconductor (S) 
can be understood in terms of normal and Andreev 
reflectioni. For systems where the effect of the Coulomb 
repulsion between electrons is negligible, the interplay 
between these two processes is governed by the trans- 
mission of the NS interface. If the transmission is high 
(close to unity) Andreev reflection dominates and the 
low-energy conductance of the interface is higher than 
the conductance measured at energies above the gap. 
On the contrary, if the transmission is low (tunneling 
regime) , Andreev reflection is strongly suppressed and so 
is the low-energy conductance of the NS interface. The 
crossover from high to low interface transparency is well 
described by the theory^ and it is in good agreement with 
experimental results on a variety of systems 3 . 

When the normal conductor is a low-dimensionality 
system, however, electron-electron interactions play an 
important role and the situation is considerably more 
comple34^S*i. In this case, electronic transport through 
the NS interface is not only determined by the trans- 
mission coefficient, but also by the strength of the inter- 
action. It is understood quite in general that electron- 
electron interactions suppress the probability for Andreev 
reflection 8 . For interacting systems, however, a simple 
and intuitive picture of the interplay between normal and 
Andreev reflection does not exist. Nevertheless, such a 
picture would be extremely useful for the qualitative in- 
terpretation of experiments. 

The purpose of this paper is to present a simple phe- 
nomenological description that captures the important 
aspects of electrical transport through a NS interface 
in the presence of repulsive electron-electron interac- 
tions in the normal conductor. To this end, we con- 
sider the simplest possible case of interacting electrons 
in a one-dimensional (ID) ballistic conductor connected 
to a superconductor with an arbitrary interface trans- 
parency. Motivated by recent theoretical studies on in- 
teracting normal mesoscopic conductors2*i2*ii, the effect 



of electron-electron interaction is included as an energy- 
dependent transmission probability T(E), which we sub- 
stitute into the theoretical expression for the current- 
voltage characteristics of a system with no electron- 
electron interaction. We substantiate the validity of the 
proposed phenomenological picture by performing ex- 
plicit calculations for the case in which T(E) is given 
by the expression valid for a barrier in an infinite ID 
interacting conductor 9 . The outcome of these calcula- 
tions reproduces qualitatively the trends found in for- 
mally correct theories^. We also illustrate how different 
aspects of our results can explain different aspects of the 
transport behavior of carbon nanotubes connected to su- 
perconducting electrodes^. 

Our work is stimulated by the description of electron 
transport through a tunnel barrier of arbitrary trans- 
mission in a ID interacting conductor, which was first 
proposed by Matveev and coworkers^. Within this pic- 
ture, the tunnel barrier induces a Friedel oscillation of 
the electron density in the one-dimensional conductor. 
In the presence of electron-electron interaction, such a 
density oscillation creates an electrostatic potential that 
contributes to scattering the electrons. Thus, to find the 
final transmission coefficient for the tunnel barrier, the 
effect of multiple reflections on the barrier and on the 
Friedel oscillation needs to be considered. It was shown 
that the net effect of the interaction is to introduce an 
energy dependence in the final transmission probability, 
such that the transmission decreases with decreasing en- 
ergy. 

This description of electron-electron interaction in 
terms of a renormalized, energy-dependent transmis- 
sion coefficient has been recently extended to deal with 
more general (i.e., not only one dimensional) mesoscopic 
conductorsiSiii. Similarly to the ID case, it was found 
that the dominant effect on low-energy electron trans- 
port is due to elastic scattering and it can be accounted 
for by an energy renormalization of the transmission co- 
efficients. At sufficiently low energy, this renormalization 
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normal metal/superconductor junction have shown that 
also in this case the effect of the interaction is to renor- 
malize the scattering coefficients^. More importantly, 
Eq. ^ captures the trend of an interaction-induced sup- 
pression of the transmission probability that quite in gen- 
eral accounts for the effect of electron-electron interaction 
on transport, as discussed above. Therefore, since our 
work only aims at providing a transparent phenomeno- 
logical description of transport across a NS interface and 
not at discussing specific details, the use of Eq.^is qual- 
itatively justified^. 

The second ingredient of our description is the expres- 
sion for the current I flowing through a NS interface of 
arbitrary transparency as a function of the applied bias 
V. This was found long ago by Blonder, Tinkham, and 
Klapwijk (BTK)£ to be: 



FIG. 1: Temperature dependence of the zero-bias differen- 
tial resistance calculated for a non-interacting (a = 0, dashed 
lines) and an interacting (a = 0.15, dotted lines) ID conduc- 
tor connected to a superconductor, for different values of the 
contact transparency (z — 0,0.22, and 1.5). For 2 = (con- 
tinuous line) the behavior is identical for the interacting and 
the non interacting case. 



results quite in general in a suppression of the transmis- 
sion coefficients with decreasing energy. These conclu- 
sions indicate that, at a qualitative level, accounting for 
the effects of electron-electron interaction in terms of a 
transmission probability decreasing with decreasing en- 
ergy is a robust result, insensitive to the details of the 
system considered. 

For the case of a tunnel barrier with (bare) trans- 
mission To hi a- n infinite ID conductor, the explicit 
interaction-induced energy dependence of T(E) reads^ 
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where < a < 1 quantifies the electron-electron inter- 
action strength (a = corresponds to no interactions) 
and Do is a high-energy cut-off determined by the en- 
ergy bandwidth of the electronic states. In what follows, 
we take T(E) as given by Eq.QJand insert it in the known 
expression for the I — V characteristics of a NS interface 
(with bare transmission To) of a non-interacting system. 

Formally this procedure is not correct, since the 
interaction-induced renormalization of the transmission 
coefficient in an infinite ID conductor differs from that 
of a tunnel barrier at a NS interface. This is because 
the presence of the superconductor changes the details of 
the Friedel oscillations and thus the specific dependence 
of T on E. Nevertheless, calculations similar to those of 
Matveev and coworkers performed for a ID interacting 



I = G Q dE[f(E-eV)-f(E)][l + A(E,z)-B(E,z)} 

J — OO 

(2) 

Here A(E, z) and B(E, z) are the energy dependent An- 
dreev and normal reflection probability, respectively, and 
Go is a constant. The parameter z is a constant that 
quantifies the amplitude of a delta-like potential barrier 
at the interface. It is related to the transmission proba- 
bility Tq as 
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FIG. 2: The dV/dl — V characteristic calculated at different 
temperatures below T c , for a high transparency contact (z = 
0.03) in the presence of interaction for the same interaction 
strength (corresponding to a = 0.6). Note that at higher 
temperature (e.g., T — 0.32 A) the dV/dl — V is similar to 
what is observed in the non-interacting case, whereas at lower 
temperature (e.g., T = 0.02 A) a zero bias peaks emerge that 
cannot be accounted for in terms of the BTK model. 
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z = corresponds to a perfectly transparent interface, 
resulting in A(E,0) = 1 and B(E,0) = for E < A 
(Andrecv limit); z>l corresponds to the tunneling limit 
in which A(E, z)«l and B(E, z) ~ 1. 

Whereas in the non-interacting case the z parameter 
is a constant, after substituting T(E) for T in Eq. [3] 
to account for the presence of interactions z acquires an 
energy dependence. We label the energy dependent z 
parameter as z e _ e . By direct substitution we have: 
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It is clear that now, even for values of T close to unity, it 
is the energy E which determines whether z e _ e is approx- 
imately or much larger than one. If, for E < A, z e _ e 
varies considerably, electron-electron interactions can in- 
duce a crossover from the Andreev limit to the tunnel- 
ing limit in an individual sample, depending on the en- 
ergy scale on which transport is probed. This crossover 
is responsible for the qualitatively different behavior of 
systems in which electron-electron interaction plays an 
important role, as compared to systems that can be de- 
scribed in terms of independent electrons. Note, in pass- 
ing, that if z = (i.e., To — 1), z e -e is also 0, implying 
that for a perfectly ballistic system, electron-electron in- 
teractions do not change the DC conductance. This is a 
well-known result that has been demonstrated from for- 
mally correct theories^ and which is reproduced by our 
phcnomenological approach. 

We illustrate the physical consequences of the energy- 
induced crossover from a transparent NS interface to the 
tunneling regime by looking at the behavior of the dif- 
ferential resistance calculated from Eq. [21 that we inte- 
grate numerically after replacing z with z e _ e . In partic- 
ular, we look at the dependence of the differential resis- 
tance on temperature and voltage for different values of 
z and a. The different dV/dl curves are normalized to 
i?jv = (1 + z 2 )/Go, the normal state resistance in the 
non-interacting case (a = 0). Note that, in general and 
contrary to the case of non-interacting electrons, the dif- 
ferential resistance of the interacting system depends on 
the bias even when eV > A, because of the energy de- 
pendence of the transmission probability. 

Fig. n shows the temperature dependence of the lin- 
ear resistance of the NS interface for different values of 
the (bare) z parameter, enabling the direct comparison 
of the non-interacting (a = 0) and the interacting case 
(with a = 0.15 in this example). As explained above, 
interactions do not have any influence if z = and the 
temperature dependence of the linear resistance is iden- 
tical in the two cases. At finite values of z, however, the 
effect of the interaction is visible. In particular, at low z 
(z = 0.22) and in the absence of interaction the resistance 
of the non-interacting system decreases with decreasing 
T below T c , so that at low temperature the resistance is 
smaller than that measured for T > T c . However, in the 
presence of interactions the behavior of the resistance is 



qualitatively different. After an initial decrease (just be- 
low T c ), the resistance increases again rapidly, so that at 
low temperature it is larger than in the normal state. 

This difference in behavior is a clear manifestation of 
the crossover from the Andreev to the tunneling limit 
mentioned above. For the interacting system, the ther- 
mal energy of the electrons at T ~ T c is sufficiently large 
to prevent a strong renormalization of the transmission 
probability at the interface, and the observed behavior 
(i.e., the decrease in resistance) is similar to the non- 
interacting case in the Andreev limit. At lower temper- 
ature, however, the renormalized transmission probabil- 
ity becomes small and the interacting system goes into 
the tunneling limit. As a consequence the resistance in- 
creases. The precise temperature at which the upturn 
occurs depends on z and on the interaction strength a, 
shifting to higher values with increasing these parame- 
ters. As shown in the figure, for sufficiently large values 
of z the resistance increases with lowering T across T c 
both in the interacting and the non-interacting case, and 
the difference in behavior becomes only of quantitative 
nature. 

A second manifestation of the crossover from Andreev 
to tunneling behavior can be seen in the dV/dl-V curves 
measured at different temperatures. Also here, the case 
exhibiting qualitative differences for the interacting and 
non-interacting cases is that of a small (bare) z parame- 
ter. The results of the calculations are shown in Fig. |3 
where z = 0.03 and a — 0.6. At higher temperature 
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FIG. 3: The dV/dl — V calculated for a fixed interface trans- 
parency [z = 0.22) and temperature T = 0.5 A, as a function 
of interaction strength (a ranging from to 0.6). Increas- 
ing the interaction results in the complete disappearance of 
the resistance suppression seen below the gap for the non- 
interacting system. 
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(approximately 0.3 A in the present case), a pronounced 
broad minimum in dV/dl is observed when the DC bias 
is smaller than the superconducting gap, as it is charac- 
teristic for Andreev reflection. As the temperature is low- 
ered, however, a sharp peak in the dV/dl curves appears 
around zero bias, causing the resistance to exceed the 
normal state value. This peak cannot be accounted for 
by the BTK-theory for non-interacting systems. Again, 
the dV/ dl peak is a consequence of the renormalization 
of the transmission probability which brings the NS in- 
terface in the tunneling regime at sufficiently low energy. 

Finally, it is instructive to look at the dependence of 
the dV/dl-V curves as a function of the electron-electron 
interaction strength a. The results of the calculations are 
shown in Fig. for the case of z — 0.22 and T = 0.5 A, 
with a ranging from to 0.6. In the non-interacting case 
(a = 0), transport occurs in the Andreev regime for all 
values of the applied bias, resulting in a suppression of 
the resistance at low energy. However, as the interaction 
strength is enhanced, a peak in dV/dl appears around 
zero bias, which becomes broader for larger values of a. 
Eventually, this peak completely dominates the behavior 
of the differential resistance for bias voltages below the 
superconducting gap. In this case, the resistance sup- 
pression characteristic of Andreev reflection is not visible 
any more. This is because, even though the bare value of 
the transmission coefficient has remained the same, the 
actual value of the transmission has changed due to the 
energy dependent renormalization induced by electron- 
electron interaction. In a non-interacting picture the ob- 
servation of such a peak in the dV/dl-V curve would be 
tentatively interpreted as due to a low transparency in- 
terface. One would however note that the shape of the 
curve is very different from what is normally observed 
in the tunneling regime, as indicated, for instance, by 
the absence of any feature at a bias corresponding to the 
superconducting gap. 

The temperature and bias voltage dependence of the 
differential resistance that we find in the interacting 
case for high (bare) values of the transmission through 



the NS interface reproduce qualitatively what has been 
found previously from formally correct calculations for 
a Superconductor-Luttinger Liquid junction!. This indi- 
cates that the phenomenological picture presented here is 
physically sound. We want to explicitly make clear, how- 
ever, that a Luttinger liquid behavior of the normal metal 
is not a necessary requirement for the observation of the 
behavior that we have described. All is needed is a suffi- 
ciently rapid suppression of the transmission probability 
with lowering energy on the scale of the superconducting 
gap. As mentioned above, in low-dimensionality samples 
such a suppression can be caused by electron-electron in- 
teraction under fairly general circumstances. 

Finally, our results exhibit qualitative agreement with 
the behavior observed experimentally in carbon nan- 
otubes contacted with superconducting electrodes. The 
case of high transmission contacts has been studied in 
Ref. [12], where an anomalous (as compared to the non 
interacting case) temperature and bias dependence of 
the differential resistance has been reported. The ex- 
perimental observations are qualitatively very similar to 
the behavior shown in Fig. ^ and |2 for the low z case. 
The case of lower transmission contacts has been studied 
in Ref. [16]. In that experiment a tunneling-like BCS- 
density of state was expected and observed without any 
clear "BCS shoulders". The absence is attributed to the 
presence of interaction which decreases the magnitude of 
this large feature. This behavior is similar to the one pre- 
dicted by our calculations, as illustrated in Fig. [3] (curves 
with a = 0.5 and 0.6). 

In conclusion, we have proposed a phenomenolog- 
ical description of the effect of electron-electron in- 
teraction on transport through a normal conduc- 
tor/superconductor interface of arbitrary transparency. 
The results obtained from this description agree qualita- 
tively with what is predicted by formally correct theories 
and reproduce the behavior observed experimentally in 
carbon nanotube/superconductor junctions. We expect 
that the simplicity of the picture proposed here will be 
useful in the interpretation of experimental data. 
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